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The conditions of quantum-classical correspondence for a system of two interacting spins are 
investigated. Differences between quantum expectation values and classical Liouville averages are 
examined for both regular and chaotic dynamics well beyond the short-time regime of narrow states. 
We find that quantum-classical differences initially grow exponentially with a characteristic exponent 
consistently larger than the largest Lyapunov exponent. We provide numerical evidence that the 
time of the break between the quantum and classical predictions scales as \og(J /K), where J is 
a characteristic system action. However, this log break-time rule applies only while the quantum- 
classical deviations are smaller than 0(h). We find that the quantum observables remain well 
approximated by classical Liouville averages over long times even for the chaotic motions of a 
few degree-of-freedom system. To obtain this correspondence it is not necessary to introduce the 
decoherence effects of a many degree-of-freedom environment. 
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I. INTRODUCTION 



' There is considerable interest in the interface between quantum and classical mechanics and the conditions that 
lead to the emergence of classical behaviour. In order to characterize these conditions, it is important to differentiate 
C^l ' two distinct regimes of quantum-classical correspondence JjJ : 

(i) Ehrenfest correspondence, in which the centroid of the wave packet approximately follows a classical trajectory. 

(ii) Liouville correspondence, in which the quantum probability distributions are in approximate agreement with those 
of an appropriately constructed classical ensemble satisfying Liouville's equation. 

Regime (i) is relevant only when the width of the quantum state is small compared to the dimensions of the system; 
if the initial state is not narrow, this regime may be absent. Regime (ii), which generally includes (i), applies to a 
much broader class of states, and this regime of correspondence may persist well after the Ehrenfest correspondence 
has broken down. The distinction between regimes (i) and (ii) has not always been made clear in the literature, 
though the conditions that delimit these two regimes, and in particular their scaling with system parameters, may be 
: S> quite different. 

The theoretical study of quantum chaos has raised the question of whether the quantum-classical break occurs 
H ' differently in chaotic states, in states of regular motion, and in mixed phase-space systems. This is well understood 
only in the case of regime (i) . There it is well-known [0-|J that the time for a minimum- uncertainty wave packet to 
expand beyond the Ehrenfest regime scales as \og{J /K) for chaotic states, and as a power of J /Ti for regular states, 
where J denotes a characteristic system action. 

The breakdown of quantum-classical correspondence, in the case of regime (ii), is less well understood, though 
it has been argued that this regime may also be delimited by a log(JY?i) break-time in classically chaotic states 
HD- Some numerical evidence in support of this conjecture has been reported in a study of the kicked rotor in the 
anomolous diffusion regime 0. (On the other hand, in the regime of quantum localization, the break-time for the 
kicked rotor seems to scale as (J /K) 2 ||.) Since the \og(J /%) time scale is rather short, it has been suggested that 
certain macroscopic objects would be predicted to exhibit non-classical behaviour on observable time scales |pj^0|- 
These results highlight the importance of investigating the characteristics of quantum-classical correspondence in more 
detail. 

In this paper we study the classical and quantum dynamics of two interacting spins. This model is convenient 
because the Hilbert space of the quantum system is finite-dimensional, and hence tractable for computations. Spin 
models have been useful in the past for exploring classical and quantum chaos P,[TT|-p"5| and our model belongs to 
a class of spin models which show promise of experimental realization in the near future p6[ . The classical limit is 
approached by taking the magnitude of both spins to be very large relative to H, while keeping their ratio fixed. For 
our model a characteristic system action is given by J ~ M, where / is a quantum number, and the classical limit is 
simply the limit of large quantum numbers, i.e. the limit I — > oo. 

In the case of the chaotic dynamics for our model, we first show that the widths of both the quantum and classical 
states grow exponentially at a rate given approximately by the largest Lyapunov exponent (until saturation at the 
system dimension). We then show that the initially small quantum-classical differences also grow at an exponential 
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rate, with an exponent X qc that is independent of the quantum numbers and at least twice as large as the largest 
Lyapunov exponent. We demonstrate how this exponential growth of differences leads to a log break-time rule, 
tb — Xqc hi(Zp/7i), delimiting the regime of Liouville correspondence. The factor p, measured in units of K, is some 
preset tolerance that defines a break between the quantum and classical expectation values. However, we also show 
that this logarithmic rule holds only if the tolerance p for quantum-classical differences is chosen extremely small, in 
particular p < 0(h). For larger values of the tolerance, the break-time does not occur on this log time-scale and may 
not occur until the recurrence time. In this sense, log break-time rules describing Liouville correspondence are not 
robust. These results demonstrate that, for chaotic states in the classical limit, quantum observables are described 
approximately by Liouville ensemble averages well beyond the Ehrenfest time-scale, after which both quantum and 
classical states have relaxed towards equilibrium distributions. This demonstration of correspondence is obtained for 
a few degree-of-frcedom quantum system of coupled spins that is described by a pure state and subject only to unitary 
evolution. 

This paper is organised as follows. In section II we describe the quantum and classical versions of our model. 
Since the model is novel we examine the behaviours of the classical dynamics in some detail. In section III we define 
the initial quantum states, which are SU(2) coherent states, and then define a corresponding classical density on 
the 2-sphere which is a good analog for these states. We show in the Appendix that a perfect match is impossible: 
no distribution on S 2 can reproduce the moments of the SU(2) coherent states exactly. In section IV we describe 
our numerical techniques. In section V we examine the quantum dynamics in regimes of classically chaotic and 
regular behaviour and demonstrate the close quantitative correspondence with the Liouville dynamics that persists 
well after the Ehrenfest break-time. In section VI we characterize the growth of quantum-classical differences in the 
time-domain. In section VII we characterize the scaling of the break-time for small quantum-classical differences and 
also examine the scaling of the maximum quantum-classical differences in the classical limit. 



II. THE MODEL 



We consider the quantum and classical dynamics generated by a non-integrable model of two interacting spins, 

oo 

H = a(S z + L z ) + cS x L x S(t - n) (1) 

n— — oo 

where S = (S x , S y , S z ) and L = (L x , L y , L z ). The first two terms in ([j]) correspond to simple rotation of both spins 
about the z-axis. The sum over coupling terms describes an infinite sequence of <5-function interactions at times t = n 
for integer n. Each interaction term corresponds to an impulsive rotation of each spin about the x-axis by an angle 
proportional to the x-component of the other spin. 



A. The Quantum Dynamics 

To obtain the quantum dynamics we interpret the Cartesian components of the spins as operators satisfying the 
usual angular momentum commutation relations, 

[Si, Sj] — iCijkSk 
[Li, Lj] — ieijkLk 

In the above we have set Ti = 1 and introduced the total angular momentum vector J = S + L. 

The Hamiltonian (Q) possesses kinematic constants of the motion, [S 2 , H] = and [L 2 , H] — 0, and the total state 
vector l^} can be represented in a finite Hilbert space of dimension (2s + I) x (21 + I). This space is spanned by the 
orthonormal vectors \s, m s ) ® \l, mi) where m s £ {s, s — I, . . . , — s} and mi £ {1,1 — 1, ... , —I}. These are the joint 
eigenvectors of the four spin operators 

S 2 \s,l,m s ,m l ) = s(s + l)\s,l,m s ,mi) 

S z \s,l,m s ,mi) = m s \s,l,m s ,mi) (2) 
L 2 |s,;,m s ,77i ; ) = 1(1 + \)\s,l,m s ,mi) 
L z \s,l,m s ,mi) = mi\s,l,m s ,mi). 
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The periodic sequence of interactions introduced by the (5-function produces a quantum mapping. The time- 
evolution for a single iteration, from just before a kick to just before the next, is produced by the unitary transfor- 
mation, 

\if>(n + l))=F\il>(n)), (3) 

where F is the single-step Floquet operator, 

F = cxp [—ia(Sz + L z )\ cxp [—icS x L x ] . (4) 

Since a is a rotation its range is 2tt radians. The quantum dynamics are thus specified by two parameters, a and c, 
and two quantum numbers, s and I. 

An explicit representation of the single-step Floquet operator can be obtained in the basis (0) by first re-expressing 
the interaction operator in (Q) in terms of rotation operators, 

exp [-icS x ® L x ] = [R (s) (9, cj>) ® R w (6, 0)] exp [-icS z <8> L„] 

x[R < - s \e,<f>)^R^(d,cl>)}- 1 , (5) 

using polar angle 9 = tt/2 and azimuthal angle = 0. Then the only non-diagonal terms arise in the expressions for 
the rotation matrices, which take the form, 

{j,m'\R^(6,4>)\j,m) = aq)(-imV)d$, m (0). (6) 

The matrix elements, 

d { 3, m (8) = (j>™'\eM-iOJy)\3,rn) (7) 

are given explicitly by Wigner's formula Jl9| . 

We are interested in studying the different time-domain characteristics of quantum observables when the cor- 
responding classical system exhibits either regular or chaotic dynamics. In order to compare quantum systems 
with different quantum numbers it is convenient to normalize subsystem observables by the subsystem magnitude 
\J (L 2 ) = yjlil + 1). We denote such normalized observables with a tilde, where 

and the normalized variance at time n is defined as, 

We are also interested in evaluating the properties of the quantum probability distributions. The probability 
distribution corresponding to the observable L z is given by the trace, 



P z {mi)=Tv\p (l >{n)\l,m l ){l,m l \\ = (Z,m,|pW (n)\l,mi), (10) 

where p^ l \n) = Tr^ [ \ijj(n))(ip(n)\ \s,m s )(s,m s \ } is the reduced state operator for the spin L at time n and Tr^ 
denotes a trace over the factor space corresponding to the spin S. 

B. Classical Map 

For the Hamiltonian (Q) the corresponding classical equations of motion are obtained by interpreting the angular 
momentum components as dynamical variables satisfying, 

{Si, Sj} = €{jkSk 

{Li, Lj} = eijkLk 
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with {■, •} denoting the Poisson bracket. The periodic 8- function in the coupling term can be used to define surfaces 
at t = n, for integer n, on which the time-evolution reduces to a stroboscopic mapping, 



nrt+1 

L" +1 



fn+l 



L 



n+l _ 



S% cos(a) - cos( 7 rL; 1 ) - S% sin( 7 rL; 1 ) sin(a), 
S 1 ™ cos( 7 rL; 1 ) - S 7 z l sin( 7 rL; l )j cos(a) + S*; 1 sin(a), 
S^co S ( 7 rL2) + S; S in( 7 rL2), 
V% cos(a) - [L; cos( 7 S£) - L n z sin( 7 S£)] sin(a), 

£J cos( 7 S£) - VI sin( 7 5J)] cos(a) + L£ sin(a), 
L"cos( 7 ^)+Z^sin( 7 ^), 



(11) 



where L = L/|L| , S = S/|S| and we have introduced the parameters 7 = c|S| and r — |L|/|S|. The mapping 
equations (11) describe the time-evolution of (0) from just before one kick to just before the next. 

Since the magnitudes of both spins are conserved, {S 2 ,H} = {h 2 ,H} = 0, the motion is actually confined to the 
four-dimensional manifold V = S 2 x S 2 , which corresponds to the surfaces of two spheres. This is manifest when the 
mapping ( O ) is expressed in terms of the four canonical coordinates x = (S z , <j> s , L z , tf>i), where 4> s = t&n(S y /S x ) and 
4>i = t&n(L y /L x ). We will refer to the mapping (^Tj) in canonical form using the shorthand notation x ,l+1 = F(x"). 

It is also useful to introduce a complete set of spherical coordinates 9 = (Q s , <fi s , 6i, <f>i) where 6 S — cos {S Z /\S\) and 
9i = cos~ 1 (L z /\L\). 

The classical flow ([ll]) on the reduced surface V still has a rather large parameter space; the dynamics are determined 
from three independent dimensionless parameters: a £ [0, 27r), 7 6 (— oo, oo), and r > 1. The first of these, a, controls 
the angle of free-field rotation about the z-axis. The parameter 7 = c|S is a dimensionless coupling strength and 
r = |L|/|S| corresponds to the relative magnitude of the two spins. 

We are particularly interested in the effect of increasing the coupling strength 7 for different fixed values of r. In 
Fig. [I] we plot the dependence of the classical behaviour on these two parameters for the case a = 5, which produces 
typical results. The data in this figure was generated by randomly sampling initial conditions on V, using the canonical 
measure, 



d/i(x) = dS z dcf) s dL z 



(12) 



and then calculating the largest Lyapunov exponent associated with each trajectory. Open circles correspond to 
regimes where at least 99% of the initial conditions were found to exhibit regular behaviour and crosses correspond 
to regimes where at least 99% of these randomly sampled initial conditions were found to exhibit chaotic behaviour. 
Circles with crosses through them (the superposition of both symbols) correspond to regimes with a mixed phase 
space. For the case a = 5 and with r held constant, the scaled coupling strength 7 plays the role of a perturbation 
parameter: the classical behaviour varies from regular, to mixed, to predominantly chaotic as | 7 | is increased from 
zero. 

The fixed points of the classical map (^) provide useful information about the parameter dependence of the classical 
behaviour and, more importantly, in the case of mixed regimes, help locate the zones of regular behaviour in the 4- 
dimensional phase space. We find it sufficient to consider only the four trivial (parameter-independent) fixed points 
which lie at the poles along the z-axis: two of these points correspond to parallel spins, (S Z ,L Z ) — ±(|S|, |L|), and 
the remaining two points correspond to anti-parallel spins, (S Z ,L Z ) = (±|S|, =p|L|). 

The stability around these fixed points can be determined from the eigenvalues of the tangent map matrix, M = 
dF/dx, where all derivatives are evaluated at the fixed point of interest. (It is easiest to derive M using the six non- 
canonical mapping equations ( |Tl| ) since the tangent map for the canonical mapping equations exhibits a coordinate 
system singularity at these fixed points.) The eigenvalues corresponding to the four trivial fixed points are obtained 
from the characteristic equation, 



- 2£ cos a + l] 2 ± f V r sin 



0, 



(13) 



with the minus (plus) sign corresponding to the parallel (anti-parallel) cases and we have suppressed the trivial factor 
(1 — £) 2 which arises since the six equations ( |ll|) are not independent. For the parallel fixed points we have the four 
eigenvalues, 



Ci.2 = cosa± - 



r 7 2 sin 2 a - 



— \j ±4 cos a 
2 



7 2 r sin 2 a 



sin 2 a(4 — 7 2 r), 
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£34 = cos a ± — \J rj 2 sin 2 a — ±4 cos a\J j 2 r sin 2 a — sin 2 a(4 — 7 2 r), (14) 



2 v ' 2 

and the eigenvalues for the anti-parallel cases, £ AP , are obtained from (]lj) through the substitution r — > — r. A fixed 
point becomes unstable if and only if |£| > 1 for at least one of the four eigenvalues. 



1, Mixed Phase Space: 7 = 1.215 



We are particularly interested in the behaviour of this model when the two spins are comparable in magnitude. 
Choosing the value r = 1.1 (with a = 5 as before), we determined by numerical evaluation that the anti-parallel 
fixed points are unstable for [--y | > 0. In the case of the parallel fixed points, all four eigenvalues remain on the unit 
circle, |£ p | = 1, for |7| < 1.42. This stability condition guarantees the presence of regular islands about the parallel 
fixed points pQ ]. In Fig. ^ we plot the trajectory corresponding to the parameters a = 5, r = 1.1, 7 = 1.215 and 
with initial condition 0(0) = (5°, 5°, 5°, 5°) which locates the trajectory near a stable fixed point of a mixed phase 
space (see Fig. [j].) This trajectory clearly exhibits a periodic pattern which we have confirmed to be regular by 
computing the associated Lyapunov exponent (Ax = 0). In contrast, the trajectory plotted in Fig. ^ is launched with 
the same parameters but with initial condition 0(0) = (20°, 40°, 160°, 130°), which is close to one of the unstable anti- 
parallel fixed points. This trajectory explores a much larger portion of the surface of the two spheres in a seemingly 
random manner. As expected, a computation of the largest associated Lyapunov exponent yields a positive number 
(X L = 0.04). 



2. Global Chaos: 7 = 2.835 



If we increase the coupling strength to the value 7 = 2.835, with a = 5 and r — 1.1 as before, then all four trivial 
fixed points become unstable. By randomly sampling V with 3 x 10 4 initial conditions we find that less than 0.1% 
of the kinematically accessible surface V is covered with regular islands (see Fig. |l|). This set of parameters produces 
a connected chaotic zone with largest Lyapunov exponent Al = 0.45. We will refer to this type of regime as one of 
'global chaos' although the reader should note that our usage of this expression differs slightly from that in 0). 



3. The Limit r » 1 



Another interesting limit of our model arises when one of the spins is much larger than the other, r > 1, We 
expect that in this limit the larger spin (L) will act as a source of essentially external 'driving' for the smaller spin 
(S). Referring to the coupling terms in the mapping (pd|), the 'driving' strength, or perturbation upon S from L, is 
determined from the product "/r = c|L|, which can be quite large, whereas the 'back-reaction' strength, or perturbation 
upon L from S, is governed only by the scaled coupling strength 7 = c|S|, which can be quite small. It is interesting 
to examine whether a dynamical regime exists where the larger system might approach regular behaviour while the 
smaller 'driven' system is still subject to chaotic motion. 

In Fig. I we plot a chaotic trajectory for r = 100 with initial condition 0(0) = (27°, 27°, 27°, 27°) which is located 
in a chaotic zone (Al = 0.026) of a mixed phase space (with a = 5 and 7 = 0.06). Although the small spin wanders 
chaotically over a large portion of its kinematically accessible shell 5 2 , the motion of the large spin remains confined 
to a 'narrow' band. Although the band is narrow relative to the large spin's length, it is not small relative to the 
smaller spin's length. The trajectories are both plotted on the unit sphere, so the effective area explored by the large 
spin (relative to the effective area covered by the small spin) scales in proportion to r 2 . 



C. The Liouville Dynamics 



We are interested in comparing the quantum dynamics generated by (g) with the corresponding Liouville dynamics 
of a classical distribution. The time-evolution of a Liouville density is generated by the partial differential equation, 

d -^ = -{p c ,H}, (15) 
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where H stands for the Hamiltonian (||) and x = (S z , (f> s ,L z ,(pi). 
The solution to (|l5| ) can be expressed in the compact form, 

Pc (x, t) = [ d^(y) 5(x - x(t, y)) p c (y, 0), (16) 

with measure d/i(y) given by ( |l2] ) and each time-dependent function x(i, y) 6 "P is solution of the equations of motion 
for (Q) with initial conditiony 6 V . This integral solution ([!(]) simply expresses that Liouville's equation (|l5|) describes 
the dynamics of a classical density p c (x, t) of points evolving in phase space under the Hamiltonian flow. We exploit 
this fact to numerically solve ( |l5|) by randomly generating initial conditions consistent with an initial phase space 
distribution p c (x, 0) and then time-evolving each of these initial conditions using the equations of motion (|ll|). We 
then calculate the ensemble averages of dynamical variables, 

(L z (n)) c = J d M (x)^p c (x,n). (17) 

by summing over this distribution of trajectories at each time step. 



D. Correspondence Between Quantum and Classical Models 



For a quantum system specified by the four numbers {a, c, s, I}, the corresponding classical parameters {a, 7, r} are 
determined if we associate the magnitudes of the classical angular momenta with the quantum spin magnitudes, 

\S\ c = y/s(s + l) 

|L| C = (18) 
This prescription produces the classical parameters, 

V s(s + 1) 

7 = cy/s(s + 1), (19) 

with a the same number for both models. 

We are interested in determining the behaviour of the quantum dynamics in the limit s — > 00 and I — > 00. This 
is accomplished by studying sequences of quantum models with s and / increasing though chosen such that the 
classical r and 7 are held fixed. Since s and / are restricted to integer (or half-integer) values, the corresponding 
classical r will actually vary slightly for each member of this sequence (although 7 can be matched exactly by 
varying the quantum parameter c). In the limit s — > 00 and I — > 00 this variation becomes increasingly small since 
r = y/l(l + l)/s(s+l) -> l/s. For convenience, the classical r corresponding to each member of the sequence of 
quantum models is identified by its value in this limit. We have examined the effect of the small variations in the 
value of r on the classical behaviour and found the variation to be negligible. 

III. INITIAL STATES 
A. Initial Quantum State 

We consider initial quantum states which are pure and separable, 

|^(0)) = |V.(0)) ®|Vi(0)>. (20) 
For the initial state of each subsystem we use one of the directed angular momentum states, 

\6,<b)=R^(0^Uj), (21) 
which correspond to states of maximum polarization in the direction (9, <fi). It has the properties: 



G 



(M|J 2 M = JCOS0 

(6,(p\J x ±iJ y \9,(j>) = je^ainO, (22) 

where j in this section refers to either I or s. 

The states (^lj) are the SU(2) coherent states, which, like their counterparts in the Euclidean phase space, are 
minimum uncertainty states pjj ; the normalized variance of the quadratic operator, 

Aj 2 = (M|J 2 M)-(MJ|M) 2 = 1 

i(i + i) C? + i)' 1 ' 

is minimised for given j and vanishes in the limit j —> oo. The coherent states \j,j) and |j, — j) also saturate the 
inequality of the uncertainty relation, 

(J 2 )(J 2 ) > (24) 
although this inequality is not saturated for coherent states polarized along other axes. 



B. Initial Classical State and Correspondence in the Macroscopic Limit 

We compare the quantum dynamics with that of a classical Liouville density which is chosen to match the initial 
probability distributions of the quantum coherent state. For quantum systems with a Euclidean phase space it 
is always possible to construct a classical density with marginal probability distributions that match exactly the 
corresponding moments of the quantum coherent state. This follows from the fact that the marginal distributions for 
a coherent state are positive definite Gaussians, and therefore all of the moments can be matched exactly by choosing 
a Gaussian classical density. For the SU(2) coherent state, however, we show in the Appendix that no classical density 
has marginal distributions that can reproduce even the low order moments of the quantum probability distributions 
(except in the limit of infinite j). Thus from the outset it is clear that any choice of initial classical state will exhibit 
residual discrepancy in matching some of the initial quantum moments. 

We have examined the initial state and dynamical quantum-classical correspondence using several different classical 
distributions. These included the vector model distribution described in the Appendix and the Gaussian distribution 
used by Fox and Elston in correspondence studies of the kicked top fl22fl . For a state polarized along the z-axis we 
chose the density, 



sin 



= C exp 
= C exp 



2sin 2 (|)) 



sm 



dJ z 



(25) 



with C = [2na 2 (l — exp(— 2er~ 2 ))] , instead of those previously considered, because it is periodic under 27r rotation. 
An initial state directed along (6 ,(j) ) is then produced by a rigid body rotation of ( p5| ) by an angle O about the 
y-axis followed by rotation with angle </> Q about the z-axis. 

The variance a 2 and the magnitude |J| C are free parameters of the classical distribution that should be chosen to 
fit the quantum probabilities as well as possible. It is shown in the Appendix that no classical density has marginal 
distributions which can match all of the quantum moments, so we concentrate only on matching the lowest order 
moments. Since the magnitude of the spin is a kinematic constant both classically and quantum mechanically, we 
choose the squared length of the classical spin to have the correct quantum value, 



\3\l = {jl) + (J 2 y) + {J 2 Z )=j{j+l). 



(26) 



For a state polarized along the z-axis, we have (J x ) — (J y ) = and (J 2 ) = (J 2 ) for both distributions as a 
consequence of the axial symmetry. Furthermore, as a consequence of (26), we will automatically satisfy the condition, 



2(J 2 ) c + (J 2 ) c =j(j + 1). 



(27) 



Therefore we only need to consider the classical moments, 
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(J z ) c = |J| G(cr 2 ) 
(J 2 ) c = |J|V G(a 2 ), 

calculated from the density ( |25| ) in terms of the remaining free parameter, a 2 , where, 

'1 + cxp(-2cr- 2 ^ 



G(a 2 ) 



1 -exp(-2cr- 2 )_ 

We would like to match both of these classical moments with the corresponding quantum values, 

(Jz) = 3, 
{Jl) = i/2, 

calculated for the coherent state (pl|). However, no choice of a 2 will satisfy both constraints. 
If we choose a 2 to satisfy (|3l|) exactly then we would obtain, 



a 2 = 



2j 



If we choose a 2 to satisfy exactly then we would obtain, 



* 2 = ^ + ^- 3 )- 



(28) 
(29) 



(30) 



(31) 
(32) 



(33) 



(34) 



(These expansions are most easily derived from the approximation G(a 2 ) ~ 1 — a 2 , which has an exponentially small 
error for large j . ) 

We have chosen to compromise between these values by fixing a 2 so that the ratio (J z )c/{J 2 )c has the correct 
quantum value. This leads to the choice, 



2V7(J+T) 2j 1.,- 



0(f 



(35) 



These unavoidable initial differences between the classical and quantum moments will vanish in the "classical" limit. 
To see this explicitly it is convenient to introduce a measure of the quantum-classical differences, 

SJ,(n) = \(J,(n))-(J z (n)) e \, (36) 

defined at time n. For an initial state polarised in direction (6, </>), the choice (p5|) produces the initial difference, 



SM0) = c ^ + o(r 2 ), 



(37) 



which vanishes as j ' ^ oo. 



IV. NUMERICAL METHODS 

We have chosen to study the time-periodic spin Hamiltonian (Q) because the time-dependence is then reduced to a 
simple mapping and the quantum state vector is confined to a finite dimensional Hilbert space. Consequently we can 
solve the exact time-evolution equations @) numerically without introducing any artificial truncation of the Hilbert 
space. The principal source of numerical inaccuracy arises from the numerical evaluation of the matrix elements of the 
rotation operator (j, m'\R(9, 4>)\j, m) — exp(—i(jmi')d^j m (6). The rotation operator is required both for calculation of 
the initial quatum coherent state, \0, (j)} — R{0, (f>)\j, m = j), and evaluation of the unitary Floquet operator. In order 
to maximise the precision of our results we calculated the matrix elements (&, (8) = (j, m'\ exp(—i8J y \j, m) using the 
recursion algorithm of Ref. ]23[ ] and then tested the accuracy of our results by introducing controlled numerical errors. 
For small quantum numbers (j < 50) we are able to confirm the correctness of our coded algorithm by comparing 
these results with those obtained by direct evaluation of Wigner's formula for the matrix elements eP, m (0). 

The time evolution of the Liouville density was simulated by numerically evaluating between 10 s and 10 9 classical 
trajectories with randomly selected initial conditions weighted according to the initial distribution (|25|). Such a large 
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number of trajectories was required in order to keep Monte Carlo errors small enough to resolve the initial normalized 
quantum-classical differences, which scale as l/8j 2 , over the range of j values we have examined. 

We identified initial conditions of the classical map as chaotic by numerically calculating the largest Lyapunov 
exponent, A^, using the formula, 

i N 

Ai = ]vE w (") (38) 

71=1 

where d(n) — |<5xj(n)| , with d(0) = 1. The differential <5x(n) is a difference vector between adjacent trajectories 
and thus evolves under the action of the tangent map <5x(n + l) = M-5x(n), where M is evaluated along some fiducial 
trajectory |2C|| . 

Since we are interested in studying quantum states, and corresponding classical distributions which have non-zero 
support on the sphere, it is also important to get an idea of the size of these regular and chaotic zones. By comparing 
the size of a given regular or chaotic zone to the variance of an initial state located within it, we can determine 
whether most of the state is contained within this zone. However, we can not perform this comparison by direct 
visual inspection since the relevant phase space is 4-dimensional. One strategy which we used to overcome this 
difficulty was to calculate the Lyapunov exponent for a large number of randomly sampled initial conditions and then 
project only those points which are regular (or chaotic) onto the plane spanned by S z = cos0 s and L z — cos If 
the variance of the initial quantum state is located within, and several times smaller than, the dimensions of a zone 
devoid of any of these points, then the state in question can be safely identified as chaotic (or regular). 



V. CHARACTERISTICS OF THE QUANTUM AND LIOUVILLE DYNAMICS 



A. Mixed Phase Space 

We consider the time-development of initial quantum coherent states ( ^T| ) evolved according to the mapping (||) 
using quantum numbers s = 140 and I = 154 and associated classical parameters 7 = 1.215, r ~ 1.1, and a = 5, 
which produce a mixed phase spacefsee Fig. |l|). The classical results are generated by evolving the the initial 
ensemble ( p5| ) using the mapping jTlf). In Fig. |E| we compare the time-dependence of the normalized quantum 
variance, AL 2 = [(L 2 ) — (L) 2 ]/l(l + 1), with its classical counterpart, AL 2 = [(L 2 ) c - (L) 2 ]/|L| 2 . Squares (diamonds) 
correspond to the dynamics of an initial quantum (classical) state centered at 9(0) — (20°, 40°, 160°, 130°), which 
is located in the connected chaotic zone near one of the unstable fixed points of the classical map. Crosses (plus 
signs) correspond to an initial quantum (classical) state centered on the initial condition 8(0) = (5°, 5°, 5°, 5°), which 
is located in the regular zone near one of the stable fixed points. For both initial conditions the quantum and 
classical results are nearly indistinguishable on the scale of the figure. In the case of the regular initial condition, the 
quantum variance remains narrow over long times and, like its classical counterpart, exhibits a regular oscillation. In 
the case of the chaotic initial condition the quantum variance also exhibits a periodic oscillation but this oscillation 
is superposed on a very rapid, approximately exponential, growth rate. This exponential growth persists until the 
variance approaches the system size, that is, when AL 2 ~ 1 . The initial exponential growth of the quantum variance 
in classically chaotic regimes has been observed previously in several models and appears to be a generic feature of 
the quantum dynamics; this behaviour of the quantum variance is mimicked very accurately by the variance of an 



initially well-matched classical distribution 1 17 22j,|24| . 

For well-localized states, in the classical case, the exponential growth of the distribution variance in chaotic zones is 
certainly related to the exponential divergence of the underlying trajectories, a property which characterizes classical 
chaos. To examine this connection we compare the observed exponential rate of growth of the widths of the classical 
(and quantum) state with the exponential rate predicted from the classical Lyapunov exponent. For the coherent 
states the initial variance can be calculated exactly, AL 2 (0) = 1/(1 + 1). Then, assuming exponential growth of this 
initial variance we get, 

AL 2 (n) ~ jexp(2X w n) for n < t saU (39) 
where a factor of 2 is included in the exponent since AL 2 corresponds to a squared length. The dotted line in Fig. 



corresponds to the prediction (39) with X w = \l = 0.04, the value of the largest classical Lyapunov exponent. As can 
be seen from the figure, the actual growth rate of the classical (and quantum) variance of the chaotic initial state is 
significantly larger than that predicted using the largest Lyapunov exponent. For comparison purposes we also plot a 
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solid line in Fig. |5| corresponding to (39) using X w = 0.13, which provides a much closer approximation to the actual 
growth rate. We find, for a variety of initial conditions in the chaotic zone of this mixed regime, that the actual 
classical (and quantum) variance growth rate is consistently larger than the simple prediction (|3^) using Al for the 
growth rate. This systematic bias requires some explanation. 

As pointed out in p2| , the presence of some discrepancy can be expected from the fact that the Lyapunov exponent 
is defined as a geometric mean of the tangent map eigenvalues sampled over the entire connected chaotic zone 
(corresponding to the infinite time limit n — *■ oo) whereas the actual growth rate of a given distribution over a small 
number of time-steps will be determined largely by a few eigenvalues of the local tangent map. In mixed regimes these 
local eigenvalues will vary considerably over the phase space manifold and the product of a few of these eigenvalues 
can be quite different from the geometric mean over the entire connected zone. 

However, we find that the actual growth rate is consistently larger than the Lyapunov exponent prediction. It is 
well known that in mixed regimes the remnant KAM tori can be 'sticky'; these sticky regions can have a significant 
decreasing effect on a calculation of the Lyapunov exponent. In order to identify an initial condition as chaotic, 
we specifically choose initial states that are concentrated away from these KAM surfaces (regular islands). Such 
initial states will then be exposed mainly to the larger local expansion rates found away from these surfaces. This 
explanation is supported by our observations that, when we choose initial conditions closer to these remnant tori, we 
find that the growth rate of the variance is significantly reduced. These variance growth rates are still slightly larger 
than the Lyapunov rate, but this is not surprising since our initial distributions are concentrated over a significant 
fraction of the phase space and the growth of the distribution is probably more sensitive to contributions from those 
trajectories subject to large eigenvalues away from the KAM boundary than those stuck near the boundary. These 
explanations are further supported by the results of the following section, where we examine a phase space regime 
that is nearly devoid of regular islands. In these regimes we find that the Lyapunov exponent serves as a much better 
approximation to the variance growth rate. 



B. Regime of Global Chaos 



If we increase the dimensionless coupling strength to 7 = 2.835, with a = 5 and r ~ 1.1 as before, then the classical 
flow is predominantly chaotic on the surface V (see Fig. [I]). Under these conditions we expect that generic initial 
classical distributions (with non-zero support) will spread to cover the full surface V and then quickly relax close to 
microcanonical equilibrium. We find that the initially localised quantum states also exhibit these generic features 
when the quantum map is governed by parameters which produce these conditions classically. 

For the non-autonomous Hamiltonian system ( pT| ) the total energy is not conserved, but the two invariants of motion 
L 2 and S 2 confine the dynamics to the 4-dimensional manifold P — S 2 x S 2 , which is the surface of two spheres. 
The corresponding microcanonical distribution is a constant on this surface, with measure (|l2]), and zero elsewhere. 
From this distribution we can calculate microcanonical equilibrium values for low order moments, where, for example, 
{L z } — (4ir)~ 2 JL L z dfj, — and {AL 2 } = {L 2 } — {L} 2 = |L| 2 . The symbols {•} denote a microcanonical average. 

To give a sense of the accuracy of the correspondence between the classical ensemble and the quantum dynamics 
in Fig. H we show a direct comparison of the dynamics of the quantum expectation value (L z ) with I = 154 and the 
classical distribution average (L z ) c for an initial coherent state and corresponding classical distribution centered at 
— (45°, 70°, 135°, 70°). To guide the eye in this figure we have drawn lines connecting the stroboscopic points of the 
mapping equations. The quantum expectation value exhibits essentially the same dynamics as the classical Liouville 
average, not only at early times, that is, in the initial Ehrenfest regime jj],^5|, but for times well into the equilibrium 
regime where the classical moment (L z ) has relaxed close to the microcanonical equilibrium value {L z } = 0. We 
have also provided results for a single trajectory launched from the same initial condition in order to emphasize the 
qualitatively distinct behaviour it exhibits. 

In Fig. ^ we show the exponential growth of the normalized quantum and classical variances on a semilog plot for 
the same set of parameters and quantum numbers. Numerical data for (a) correspond to initial condition 0(0) = 
(20°, 40°, 160°, 130°) and those for (b) correspond to 0(0) = (45°, 70°, 135°, 70°). As in the mixed regime case, the 
quantum-classical differences are nearly imperceptible on the scale of the figure, and the differences between the 
quantum and classical variance growth rates are many orders of magnitude smaller than the small differences in the 
growth rate arising from the different initial conditions. 

In contrast with the mixed regime case, in this regime of global chaos the prediction (^9|) with A u , = Xl = 0.45 
now serves as a much better approximation of the exponential growth rate of the quantum variance, and associated 
relaxation rate of the quantum and classical states. In this regime the exponent A^, is also much larger than in the 
mixed regime case due to the stronger degree of classical chaos. As a result, the initially localised quantum and 
classical distributions saturate at system size much sooner. 
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It is useful to apply ( |39} ) to estimate the time-scale at which the quantum (and classical) distributions saturate at 
system size. From the condition AL 2 (i sat ) ~ 1 and using (139) we obtain, 



t sat ^(2X w )- 1 \n(l) (40) 

which serves as an estimate of this characteristic time-scale. In the regimes for which the full surface V is predominately 
chaotic, we find that the actual exponential growth rate of the width of the quantum state, X w , is well approximated 
by the largest Lyapunov exponent A^. For a = 5 and r = 1.1, the approximation X w ~ Al holds for coupling strengths 
7 > 2, for which more than 99% of the surface V is covered by one connected chaotic zone (see Fig. |l|). 

By comparing the quantum probability distribution to its classical counterpart, we can learn much more about 
the relaxation properties of the quantum dynamics. In order to compare each mi value of the quantum distribution, 
Pz(jni), with a corresponding piece of the continuous classical marginal probability distribution, 



p c {l z ) = ds z d4> s d4>i Pc (e s ,(j> s A^i), (4i) 

we discretize the latter into 2j + 1 bins of width % = 1. This procedure produces a discrete classical probability 
distribution P°(mj) which prescribes the probability of finding the spin component L z in the interval [mj+1/2, mi — 1/2] 
along the z-axis. 

To illustrate the time-development of these distributions we compare the quantum and classical probability distri- 
butions for three successive values of the kick number n, using the same quantum numbers and initial condition as in 
Fig. [| In Fig. || the initial quantum and classical states are both well-localised and nearly indistinguishable on the 
scale of the figure. At time n — 6 ~ t sa t, shown in Fig. ^, both distributions have grown to fill the accessible phase 
space. It is at this time that the most significant quantum-classical discrepancies appear. 

For times greater than t sat , however, these emergent quantum-classical discrepencies do not continue to grow, since 
both distributions begin relaxing towards equilibrium distributions. Since the dynamics are confined to a compact 
phase space, and in this parameter regime the remnant KAM tori fill a negligibly small fraction of the kinematicaly 
accessible phase space, we might expect the classical equilibrium distribution to be very close to the microcanonical 
distribution. Indeed such relaxation close to microcanonical equilibrium is apparent for both the quantum and the 
classical distribution at very early times, as demonstrated in Fig. p0[, corresponding to n = 15. 

Thus the signature of a classically hyperbolic flow, namely, the exponential relaxation of an arbitrary distribution 
(with non-zero measure) to microcanonical equilibrium [ p6[ , holds to good approximation in this model in a regime of 
global chaos. More suprisingly, this classical signature is manifest also in the dynamics of the quantum distribution. In 
the quantum case, however, as can be seen in Fig. [l0|, the probability distribution is subject to small irreducible time- 
dependent fluctuations about the classical equilibrium. We examine these quantum fluctuations in detail elsewhere 



VI. TIME-DOMAIN CHARACTERISTICS OF QUANTUM-CLASSICAL DIFFERENCES 

We consider the time dependence of quantum-classical differences defined along the z-axis of the spin L, 

5L z {n) = \{L z {n))-{L z (n)) e \, (42) 

at the stroboscopic times t = n. In Fig. |Il]we compare the time-dependence of SL z (n) on a semi-log plot for a chaotic 
state (filled circles), with 0(0) = (20°, 40°, 160°, 130°), and a regular state (open circles), 0(0) = (5°, 5°, 5°, 5°), evolved 
using the same mixed regime parameters (7 = 1.215 and r ~ 1.1) and quantum numbers (I = 154) as in Fig. 

We are interested in the behaviour of the upper envelope of the data in Fig. |ll]. For the regular case, the upper 
envelope of the quantum-classical differences grows very slowly, as some polynomial function of time. For the chaotic 
case, on the other hand, at early times the difference measure (^2j) grows exponentially until saturation around n = 15, 
which is well before reaching system dimension, |L| ~ I = 154. After this time, which we denote t*, the quantum- 
classical differences exhibit no definite growth, and fluctuate about the equilibrium value SL Z ~ 1 <C |L|. In Fig. |ll| 
we also include data for the time-dependence of the Ehrenfest difference \(L Z ) — L z \, which is defined as the difference 
between the quantum expectation value and the dynamical variable of a single trajectory initially centered on the 
quantum state. In contrast to 6L Z , the rapid growth of the Ehrenfest difference continues until saturation at the 
system dimension. 

In Fig. |l2] we compare the time-dependence of the quantum-classical differences in the case of the chaotic initial 
condition 0(0) = (20°, 40°, 160°, 130°) for quantum numbers I — 22 (filled circles) and I = 220 (open circles), using 
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the same parameters as in Fig. |Tl|. This demonstrates the remarkable fact that the exponential growth terminates 
when the difference measure reaches an essentially fixed magnitude (SL Z ~ 1 as for the case I — 154), although the 
system dimension differs by an order of magnitude in the two cases. 

In Fig. [l3| we consider the growth of the quantum-classical difference measure SL z (n) in a regime of global chaos, 
for I = 154, and using the same set of parameters as those examined in Fig. M (7 = 2.835 and r ~ 1.1). Again the 
upper envelope of the difference measure SL z (n) exhibits exponential growth at early times, though in this regime of 
global chaos the exponential growth persists only for a very short duration before saturation at t* ~ 6. The initial 
condition 9(0) = (20°, 40°, 160°, 130°) is a typical case (filled circles), where, as seen for the mixed regime parameters, 
the magnitude of the difference at the end of the exponential growth phase saturates at the value SL z (t*) ~ 1, which 
does not scale with the system dimension (see Fig. |lf|). The initial condition 9(0) — (45°, 70°, 135°, 70°) (open circles) 
leads to an anomolously large deviation at the end of the exponential growth phase, SL z (t*) ~ 10, though still small 
relative to the system dimension |L| ~ 154. This deviation is transient however, and at later times the magnitude of 
quantum-classical differences fluctuates about the equilibrium value 6L Z ~ 1. The quantum-classical differences are a 
factor of l/l smaller than typical differences between the quantum expectation value and the single trajectory, which 
are of order system dimension (see Fig. |6|) as in the mixed regime case. 

In all cases where the initial quantum and classical states are launched from a chaotic zone we find that the initial 
time-dependence of quantum-classical differences compares favorably with the exponential growth ansatz, 

5L z (n) ~ -- exp(X qc n) for n<t*, (43) 
of 

where the exponent X qc is a new exponent subject to numerical measurement Jl7| . The prefactor 1/81 is obtained from 
( |37| ) though we have dropped the cos 9 factor that specifies the exact initial difference for L z . Since contributions from 
the initial differences in other mismatched moments will generally mix under the dynamical flow, it is appropriate to 
consider an effective initial difference for the prefactor in (^3|). The prefactor 1/81 is obtained by accounting for the 
initial contributions from the 3 cartesian components, [S 2 L x (0) + S 2 L y (0) + 5 2 L Z (0] 1 ^ 2 = 1/81. 

We are interested in whether the Lyapunov exponent A^ is a good approximation to X qc . In Fig. [l^ we plot ( |4^ ) 
with X qc = Xl = 0.04 (dotted line) for I = 220. Clearly the largest Lyapunov exponent severly underestimates 
the exponential growth rate of the quantum-classical differences, in this case by more than an order of magnitude. 
The growth rate of the state width, X w = 0.13 , is also several times smaller than the initial growth rate of the 
quantum-classical differences. In the case of Fig. [l3[ corresponding to a regime of global chaos with a much larger 
Lyapunov exponent, we plot ( |43| ) with X qc = Xl = 0.45 (dotted line), demonstrating that, in this regime too the 
largest Lyapunov exponent underestimates the initial growth rate of the quantum-classical difference measure 5L Z (n) . 



We also find, from inspection of our results, that the time t* at which the exponential growth (43) terminates can 
be estimated from t sat , the time-scale on which the distributions saturate at or near system size (40|). In the case 



of the chaotic initial condition of Fig 
This should be compared with Fig 



for which 7 = 1.215, visual inspection of the figure suggests that t sat ~ 18. 
where the exponential growth of 6L z (n) ends rather abruptly at t* ~ 15. In 
Fig. fj], corresponding to a regime of global chaos (7 = 2.835), the variance growth saturates much earlier, around 
t sat ~ 6 for both initial conditions. From Fig. [l^ we can estimate that the initial exponential growth of the quantum- 
classical differences for these two initial conditions also ends around n ~ 6. As we increase 7 further, we find that the 
exponential growth phase of quantum-classical differences SL z (n) is shortened, lasting only until the corresponding 
quantum and classical distributions saturate at system size. For 7 ~ 12, with Xl — 1.65, the chaos is sufficiently 
strong that the initial coherent state for I = 154 spreads to cover V within a single time-step. Similarly the initial 
difference measure 8L Z (0) ~ 0.001 grows to the magnitude 8L Z (1) ~ 1 within a single time-step and subsequently 
fluctuates about that equilibrium value. We have also inspected the variation of t* with the quantum numbers and 
found it to be consistent with the logarithmic dependence of t sat in (Efl) . 



VII. CORRESPONDENCE SCALING IN THE CLASSICAL LIMIT 

We have assumed in ( [43| ) that the exponent X qc is independent of the quantum numbers. A convenient way of 
confirming this, and also estimating the numerical value of X qc , is by means of a break-time measure. The break-time 
is the time tb(l,p) at which quantum-classical differences exceed some fixed tolerance p, with the classical parameters 
and initial condition held fixed. Setting 5L z (tb) = p in (Q), we obtain % in terms of p, I and X qc , 

t b ~ A~ c x ln(8 p I) provided p < 0(1). (44) 

The restriction p < 0(1), which plays a crucial role in limiting the robustness of the break-time measure (|44|), is 
explained and motivated further below. 
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The explicit form we have obtained for the argument of the logarithm in (|44|) is a direct result of our estimate 
that the initial quantum-classical differences arising from the Cartesian components of the spin provide the dominant 
contribution to the prefactor of the exponential growth ansatz ([h}]). Differences in the mismatched higher order 
moments, as well as intrinsic differences between the quantum dynamics and classical dynamics, may also contribute 
to this effective prefactor. Wc have checked that the initial value SL z (0) ~ 1/8/ is an adequate estimate by comparing 
the intercept of the quantum-classical data on a semilog plot with the prefactor of ( fl3| ) for a variety of I values (see 
e.g. Fig. fl 

In Fig. [14| we examine the scaling of the break-time for I values ranging from 11 to 220 and with fixed tolerance 
p = 0.1. The break-time can assume only the integer values t = n and thus the data exhibits a step- wise behaviour. For 
the mixed regime parameters, 7 = 1.215 and r ~ 1.1 (filled circles), with initial condition 6(0) = (20°, 40°, 160°, 130°), 
a non-linear least squares fit to @ gives X qc = 0.43. This fit result is plotted in the figure as a solid line. The 
close agreement between the data and the fit provides good evidence that the quantum-classical exponent X qc is 
independent of the quantum numbers. To check this result against the time-dependent 5L z (n) data, we have plotted 
the exponential curve © with X qc = 0.43 in Fig. [ll] using a solid line and in Fig. [ll] using a solid line for I = 22 and 
a dotted line for / = 220. The exponent obtained from fitting ( |44] ) serves as an excellent approximation to the initial 
exponential growth ([h|) of the quantum-classical differences in each case. 

In Fig. 14 we also plot break-time results for the global chaos case 7 = 2.835 and r ~ 1.1 (open circles) with initial 
condition 6(0) — (45°, 70°, 135°, 70°). In this regime the quantum-classical differences grow much more rapidly and, 
consequently, the break-time is very short and remains nearly constant over this range of computationally accessible 
quantum numbers. Due to this limited variation, in this regime we can not confirm (|44|), although the data is consistent 
with the predicted logarithmic dependence on /. Moreover, the break-time results provide an effective method for 
estimating X qc if we assume that (^4|) holds. The same fit procedure as detailed above yields the quantum-classical 
exponent X qc = 1.1. This fit result is plotted in Fig. [l4]as a solid line. More importantly, the exponential curve J43|), 
plotted with fit result X qc = 1.1, can be seen to provide very good agreement with the initial growth rate of Fig. |13| 
for either initial condition, as expected. 

In the mixed regime (7 = 1.215), the quantum-classical exponent X qc = 0.43 is an order of magnitude greater than 
the largest Lyapunov exponent Xl = 0.04 and about three times larger than the growth rate of the width X w = 0.13. 
In the regime of global chaos (7 = 2.835) the quantum-classical exponent X qc = 1.1 is a little more than twice as large 
as the largest Lyapunov exponent Xl = 0.45. 

The condition p < 0(1) is a very restrictive limitation on the domain of application of the log break-time ( ffi] ) and 
it is worthwhile to explain the significance of this restriction. In the mixed regime case of Fig. [ll], with I = 154, we 
have plotted the tolerance values p = 0.1 (dotted line) and p = 15.4 (sparse dotted line). The tolerance p = 0.1 is 
exceeded at t = 11, while the quantum-classical differences are still growing exponentially, leading to a log break-time 
for this tolerance value. For the tolerance p = 15.4 <C |L|, on the other hand, the break-time does not occur on a 
measurable time-scale, whereas according to the logarithmic rule (Elf), with I = 154 and X qc = 0.43, we should expect 
a rather short break-time if, ~ 23. Consequently the break-time (|44|), applied to delimiting the end of the Liouville 
regime, is not a robust measure of quantum-classical correspondence. 

Our definition of the break-time ( f44| ) requires holding the tolerance p fixed in absolute terms (and not as fraction 
of system dimension as in ||) when comparing systems with different quantum numbers. Had we chosen to compare 
systems using a fixed relative tolerance, /, then the break-time would be of the form if, ~ X~^ ln(8 / I 2 ) and subject 
to the restriction / < 0(1/ 1). Since / — > in the classical limit, this form emphasizes that the log break-time applies 
only to differences that are vanishing fraction of the system dimension in that limit. 

Although we have provided numerical evidence (in Fig. [l2] ) of one mixed regime case in which the largest quantum- 
classical differences occuring at the end of the exponential growth period remain essentially constant for varying 
quantum numbers, SL z (t*) ~ 0(1), we find that this behaviour represents the typical case for all parameters and 
initial conditions which produce chaos classically. To demonstrate this behaviour we consider the the scaling (with 
increasing quantum numbers) of the maximum values attained by SL z (n) over the first 200 kicks, 8L™ ax . Since 
t* <C 200 over the range of I values examined, the quantity i5L™ aa: is a rigorous upper bound for 5L z (t*). 

In Fig. [[5] we compare SL™ ax for the two initial conditions of Fig. jl^ and using the global chaos parameters 



(7 = 2.835, r ~ 1.1). The filled circles in Fig. |15| correspond to the initial condition 6(0) = (20°, 40°, 160°, 130°). As 
in the mixed regime, the maximum deviations exhibit little or no scaling with increasing quantum number. This is 
the typical behaviour that we have observed for a variety of different initial conditions and parameter values. These 
results motivate the generic rule, 

SL z (t*) < SL r z nax ~ 0(1/0- (45) 

Thus the magnitude of quantum-classical differences reached at the end of the exponential growth regime, expressed 
as a fraction of the system dimension, approaches zero in the classical limit. 
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However, for a few combinations of parameters and initial conditions we do observe a 'transient' discrepancy peak 
occuring at t ~ t* that exceeds 0(1). This peak is quickly smoothed away by the subsequent relaxation of the quantum 
and classical distributions. This peak is apparent in Fig. |l3| (open circles), corresponding to the most conspicuous case 
that we have identified. This case is apparent as a small deviation in the normalized data of Fig. ||. The scaling of the 



magnitude of this peak with increasing I is plotted with open circles in Fig. 15, The magnitude of the peak initially 
increases rapidly but appears to become asymptotically independent of I. The other case that we have observed occurs 
for the classical parameters 7 = 2.025, with r ~ 1.1 and a = 5, and with initial condition 8(0) = (20°, 40°, 160°, 130°). 
We do not understand the mechanism leading to such transient peaks, although they are of considerable interest since 
they provide the most prominent examples of quantum-classical discrepancy that we have observed. 



VIII. DISCUSSION 



In this study of a non-integrable model of two interacting spins we have characterized the correspondence between 
quantum expectation values and classical ensemble averages for intially localised states. We have demonstrated that in 
chaotic states the quantum-classical differences initially grow exponentially with an exponent \ qc that is consistently 
larger than the largest Lyapunov exponent. In a study of the moments of the Hcnon-Heiles system, Ballcntine and 



McRae |17 lSj] have also shown that quantum-classical differences in chaotic states grow at an exponential rate with 
an exponent larger than the largest Lyapunov exponent. This exponential behaviour appears to be a generic feature 
of the short-time dynamics of quantum-classical differences in chaotic states. 

Since we have studied a spin system, we have been able to solve the quantum problem without truncation of the 
Hilbert space, subject only to numerical roundoff, and thus we are able to observe the dynamics of the quantum- 
classical differences well beyond the Ehrenfest regime. We have shown that the exponential growth phase of the 
quantum-classical differences terminates well before these differences have reached system dimension. We find that 
the time-scale at which this occurs can be estimated from the time-scale at which the distribution widths approach 
the system dimension, t sat ~ (2AU,) -1 m(Z) for initial minimum uncertainty states. Due to the close correspondence 
in the growth rates of the quantum and classical distributions, this time-scale can be estimated from the classical 
physics alone. This is useful because the computational complexity of the problem does not grow with the system 
action in the classical case. Moreover, we find that the exponent X w can be approximated by the largest Lyapunov 
exponent when the kinematic surface is predominantly chaotic. 

We have demonstrated that the exponent \ qc governing the initial growth rate of quantum-classical differences 
is independent of the quantum numbers, and that the effective prefactor to this exponential growth decreases as 
1/7. These results imply that a log break-time rule (44) delimits the dynamical regime of Liouville correspondence. 
However, the exponential growth of quantum-classical differences persists only for short times and small differences, 
and thus this log break-time rule applies only in a similarly restricted domain. In particular, we have found that the 
magnitude of the differences occuring at the end of the initial exponential growth phase does not scale with the system 
dimension. A typical magnitude for these differences, relative to the system dimension, is 0(1/1). Therefore, log(Z) 
break-time rules characterizing the end of the Liouville regime are not robust, since they apply to quantum-classical 
differences only in a restricted domain, i.e. to relative differences that are smaller than 0(1/ 1). 

This restricted domain effect does not arise for the better known log break-time rules describing the end of the 
Ehrenfest regime |l| |j. The Ehrenfest log break-time remains robust for arbitrarily large tolerances since the corre- 
sponding differences grow roughly exponentially until saturation at the system dimension p2| , p4{ . Consequently, a 
log(Z) break-time indeed implies a breakdown of Ehrenfest correspondence. However, the logarithmic break-time rule 
characterizing the end of the Liouville regime does not imply a breakdown of Liouville correspondence because it 
does not apply to the observation of quantum-classical discrepancies larger than 0(1/1). The appearance of residual 
0(1/1) quantum-classical discrepancies in the description of a macroscopic body is, of course, consistent with quantum 
mechanics having a proper classical limit . 

We have found, however, that for certain exceptional combinations of parameters and initial conditions there are 
relative quantum-classical differences occuring at the end of the exponential growth phase that can be larger than 
0(1/1), though still much smaller than the system dimension. In absolute terms, these transient peaks seem to 
grow with the system dimension for small quantum numbers but become asymptotically independent of the system 
dimension for larger quantum numbers. Therefore, even in these least favorable cases, the fractional differences 
between quantum and classical dynamics approach zero in the limit I — > 00. This vanishing of fractional differences 
is sufficient to ensure a classical limit for our model. 

Finally, contrary to the results found in the present model, it has been suggested that a log break-time delimiting 
the Liouville regime implies that certain isolated macroscopic bodies in chaotic motion should exhibit non-classical 
behaviour on observable time scales. However, since such non-classical behaviour is not observed in the chaotic motion 
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of macroscopic bodies, it is argued that the observed classical behaviour emerges from quantum mechanics only when 
the quantum description is expanded to include interactions with the many degrees-of-freedom of the ubiquitous 
environment [p|JlC|]. (This effect, called decoherence, rapidly evolves a pure system state into a mixture that is 
essentially devoid of non-classical properties.) However, in our model classical behaviour emerges in the macroscopic 
limit of a simple few degree-of-freedom quantum system that is described by a pure state and subject only to unitary 
evolution. Quantum-classical correspondence at both early and late times arises in spite of the log break-time because 
this break-time rule applies only when the quantum-classical difference threshold is chosen smaller than 0{h). In this 
sense we find that the decoherence effects of the environment are not necessary for correspondence in the macroscopic 
limit. Of course the effect of decoherence may be experimentally significant in the quantum and mesoscopic domains, 
but it is not required as a matter of principle to ensure a classical limit. 
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X. APPENDIX 



Ideally we would like to construct an initial classical density that reproduces all of the moments of the initial 
quantum coherent states. This is possible in a Euclidean phase space, in which case all Weyl-ordered moments of 
the coherent state can be matched exactly by the moments of a Gaussian classical distribution. However, below we 
prove that no classical density p c (9,(f>) that describes an ensemble of spins of fixed length |J| can be constructed 
with marginal distributions that match those of the SU(2) coherent states (pi]). Specifically, we consider the set of 
distributions on S 2 with continuous independent variables 9 £ [0,7r] and <fi £ [0, 2w), measure dp = sin0d9d(f), and 
subject to the usual normalization, 

/ dfxp c (d,cf>) = l. (46) 
Js 2 

For convenience we choose the coherent state to be polarized along the positive z-axis, p — This state is 

axially symmetric: rotations about the z-axis by an arbitrary angle <f> leave the state operator invariant. Consequently 
we require axially symmetry of the corresponding classical distribution, 

p c (6,<p)=p c (6). (47) 

We use the expectation of the quadratic operator, (J 2 ) = j(j + 1), to fix the length of the classical spins, 



I J| = VWh = (48) 

Furthermore, the coherent state \j,j) is an eigenstate of J z with moments along the z-axis given by (J™) = j n for 
integer n. Therefore we require that the classical distribution produces the moments, 

(J?) c =i". (49) 

These requirements are satisfied by the (5-function distribution, 

Pv{V) = ~~z ;— 7-, (50) 

2n sin o 

where cos6* = j/|J| defines 8 . This distribution is the familiar vector model of the old quantum theory corresponding 
to the intersection of a cone with the surface of the sphere. 

However, in order to derive an inconsistency between the quantum and classical moments we do not need to assume 
that the classical distribution is given explicitly by (|50|); we only need to make use of the the azimuthal invariance 
condition (fl7|), the length condition (^8|), and the first two even moments of (^9|). 

First we calculate some of the quantum coherent state moments along the x-axis (or any axis orthogonal to z), 
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(J™) = for odd m 

(^> = 3j 2 /4-j/4. 
In the classical case, these moments are of the form, 

(J?) c = [dJ z [ dct>p c (9)\J\ m cos m (<£)sm ro (0). (51) 



For m odd the in teg ral over <j> vanishes, as required for correspondence with the odd quantum moments. For m even 
we can evaluate (^l|) by expressing the r.h.s. as a linear combination of the z-axis moments ([l9| ) of equal and lower 
order. For m = 2 this requires substituting sin 2 (#) = 1 — cos 2 (#) into (|5l|) and then integrating over <f> to obtain 



(J 2 ) c = tt / dJ zPc (6)\J\ 2 -ir J dJ zPc (9)\3\ 2 cos 2 (9) 
= \J\/2-(J 2 )/2. 

Since (J 2 ) is determined by ( f49| ) and the length is fixed from ([l8]) we can deduce the classical value without knowing 
P(0), 

{Jl)c = (52) 
This agrees with the value of corresponding quantum moment. For m = 4, however, by a similar procedure we deduce 

(Jt) c = 3j 2 /8, (53) 
that differs from the quantum moment (J*) by the factor, 

6J* = \(J*)-{J±) c \ = \3f/8-j/4\, (54) 
concluding our proof that no classical distribution on S 2 can reproduce the quantum moments. 
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FIG. 1. Behaviour of the classical mapping for different values of r = |L|/|S| and 7 = c|S| with a = 5. Circles correspond to 
parameter values for which at least 99% of the surface area V produces regular dynamics and crosses correspond to parameter 
values for which the dynamics are at least 99% chaotic. Superpositions of circles and crosses correspond to parameter values 
which produce a mixed phase space. We investigate quantum-classical correspondence for the parameter values 7 = 1.215 
(mixed regime) and 7 = 2.835 (global chaos), with r = 1.1, which are indicated by filled circles. 





FIG. 2. Stroboscopic trajectories on the unit sphere launched from a regular zone of the mixed regime with 7 = 1.215, 
r = 1.1, a = 5 and (9(0) = (5°, 5°, 5°, 5°). 
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FIG. 3. Same parameters as Fig. g but the trajectory is launched from a chaotic zone of the mixed regime with initial 
condition 0(0) = (20°, 40°, 160°, 130°). 




FIG. 4. A chaotic trajectory for mixed regime parameters 7 = 0.06, r = 100, and a = 5 with 0(0) = (27°, 27°, 27°, 27°). The 
motion of the larger spin appears to remain confined to a narrow band on the surface of the sphere. 
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FIG. 5. Growth of normalized quantum and classical variances in a chaotic zone (a) and a regular zone (b) of the mixed 
phase space regime 7 = 1.215 and r ~ 1.1 with I = 154. Quantum and and classical results are nearly indistinguishable on 
this scale. In the chaotic case, the approximately exponential growth of both variances is governed by a much larger rate, 
A„ar = 0.13 (solid line), than that predicted from the largest Lyapunov exponent, \l = 0.04 (dotted line). 
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FIG. 6. Comparison of quantum expectation value and corresponding classical average (L z } c in the regime of global chaos 
7 = 2.835 and r ~ 1.1 with I = 154 and initial condition o = (45°, 70°, 135°, 70°). The points of the stroboscopic map are 
connected with lines to guide the eye. The quantum expectation value and the Liouville average exhibit esentially the same 
rate of relaxation to microcanonical equilibrium, a behaviour which is qualitatively distinct from that of the single trajectory. 
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FIG. 7. Growth of normalized quantum and classical variances in the regime of global chaos, 7 = 2.835 and r ~ 1.1 with 
1=154, for the two initial conditions cited in the text. Quantum-classical differences are nearly imperceptible on this scale. In 
this regime the largest Lyapunov exponent \l = 0.45 provides a much better estimate of the initial variance growth rate. 
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FIG. 8. Initial probability distributions for L z for 0(0) = (45°, 70°, 135°, 70°) with I = 154. The quantum and classical 
distributions are initially indistinguishable on the scale of the figure. 
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FIG. 9. Same as Fig. |8| but the states have evolved to n = 6 in the regime of global chaos 7 = 2.835 and r ~ 1.1. Both 
the quantum and classical distribution have spread to system dimension and exhibit their largest differences on this saturation 
time-scale. 
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FIG. 10. Same as Fig. but for n = 15. Both quantum and classical distributions have relaxed close to the microcanonical 
equilibrium. 
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FIG. 11. Time-dependence of quantum-classical differences in a regular zone (open circles) and a chaotic zone (filled circles) 
of mixed regime (7 = 1.215 and r ~ 1.1) with I — 154. For the chaotic state SL Z — \{L Z ) — (L z ) c \ is compared with the 
Ehrenfest difference, \ {L Z ) —L g \, between the quantum expectation value and a single trajectory (plus signs), which grows until 
saturation at system dimension. The solid line corresponds to (^) using X qc = 0.43. The horizontal lines indicate two different 
values of the difference tolerance p which may be used to determine the break-time; for p = 0.1 (dotted line) tt, occurs on a 
logarithmic time-scale, but for p — 15.4 (sparse dotted line) it, is not defined over numerically accessible time-scales. 
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FIG. 12. Growth of the quantum-classical difference 8L Z in the chaotic zone of a mixed re gim e, 7 = 1.215 and r ~ 1.1, with 
I — 22 (filled circles) and I = 220 (open circles). For I = 220 the exponential growth rate (J43) is plotted using the classical 
Lyapunov exponent, \l = 0.04 (sparse dotted line), and for both I values (^3|) is plotted using the exponent X qc = 0.43 (solid 
line for I — 22, dotted line for / = 220), which is obtained from a fit of H44) to the corresponding break-time data in Fig. [141 
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FIG. 13. Growth of quantum-classical differences in the regime of global chaos 7 = 2.835 and r ~ 1.1 with / = 154 for 
the two initial conditions cited in text. The exponential growth rate (E3J) is plotted using the classical Lyapunov exponent, 
Xl = 0.45 (dotted line), and the exponent X qc = 1.1 (solid line), which is obtained from a fit of ( p4| ) to the corresponding 
break-time data in Fig. JL4L 




FIG. 14. Scaling of the break-time using tolerance p — 0.1 as a function of increasing quantum number for the mixed 
regime parameters 7 = 1.215 and r ~ 1.1 with (9(0) = (20°, 40°, 160°, 130°) (filled circles) and for the global chaos parameters 
7 = 2.835 and r ~ 1.1 with 6(0) = (45°, 70°, 135°, 70°) (open circles). We also plot the results of fits to the log rule @, which 
produced exponents A 9C = 0.43 for 7 = 1.215 and X qc = 1.1 for 7 = 2.835. 



24 



100 



o ° o o o 



10 r 



5U: 



OOO 



0.1 



50 100 150 200 250 



FIG. 15. Maximum quantum-classical difference occuring over the first 200 kicks in the regime of global chaos (7 = 2.835, 
r ~ 1.1) plotted against increasing quantum number. These maximum values provide an upper bound on 5L z (t*) for each I. 
The data corresponding to the initial condition #(0) = (20°, 40°, 160°, 130°) (filled circles) represent a typical case in which the 
maximum quantum-classical differences do not vary significantly with I. The large deviations observed for the initial condition 
8(0) = (45°, 70°, 135°, 70°) (open circles) are an exceptional case, with maximum differences growing rapidly for small quantum 
numbers but tending asymptotically toward independence of I. These curves provide an upper bound on the tolerance values 
p for which the break-time measure scales logarithmicly with I. 
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